Charge ordering in quasi-one-dimensional systems with frustrating interactions 
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Motivated by the co-existing charge and spin order found in strongly correlated ladder systems, 
we study an effective pseudospin model on a coupled two-leg ladder. A bosonisation analysis yields 
a rich phase diagram showing Wigner/Peierls charge order and Neel/dimer spin order. In a broad 
parameter regime, the orbital antiferromagnetic phase is found to be stable. An intermediate gapless 
phase of finite width is found to lie in between two charge-ordered gapped phases. Our work is 
potentially relevant for a unified description of a class of strongly correlated, quarter-filled ladder 
systems. 
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Strongly correlated ladder systems are fascinating can- 
didates for studying the interplay of spin and charge 
ordering, and their combined influence on the emer- 
gence of novel ground states 0. Several recent stud- 
ies provide experimental realisations of such systems, 
exhibiting diverse phenomena like charge order (CO), 
antiferromagnetism (AF) and unconventional supercon- 
ductivity (uSC) as functions of suitable control param- 
eters |2j, |3|. The existence of several non-perturbative 
theoretical techniques in one dimension have also resulted 
in the study of the emergence of exotic phases from in- 
stabilities of the high-T Luttinger liquid Q . Given that 
these systems are Mott insulators, longer-range Coulomb 
interactions are relevant in understanding CO/AF/uSC 
phases. Inspite of some recent work a detailed un- 
derstanding of the effects of longer-range interactions in 
quasi-lD systems is a largely unexplored problem. Fur- 
ther, attention has mostly focussed on studies of models 
at 1/2-nllmg 0. 

Here, we study an effective pseudospin model describ- 
ing charge degrees of freedom of a one-dimensional 1/4- 
filled electronic system. Such an effective model can be 
derived from an extended Hubbard model (with longer- 
range interactions) in a number of physically relevant 
cases 0, 0, ll| ■ It has also been suggested that this 
may be a relevant starting point for the study of charge- 
order in organics. While the weak coupling limit of the 
underlying extended Hubbard model has been recently 
studied [10, lUJ| , we note that the real systems under con- 
sideration are generically in the strong coupling regime 
of the model [aHOl- To the best of our knowledge, this 
regime has not been studied in sufficient detail. 

We start with the effective pseudospin Hamiltonian 



(1) 



where 2t is the transverse field, V is the n.n interaction 
pseudospin coupling. Here, the r represent charge de- 
grees of freedom in an effective model derived from a 
more basic electronic Hamiltonian for a 1 / 4-filled system 
This is the ID Ising model in a transverse field. 



Rotating the pseudospin axis t x — > t z , r z — > ~t x and 
introducing a bond-fermion repulsion U± J2i a b=ia( ni , a ~ 
rii+\ t a)(ni,b — rii+i.f)) as well as a bond-fermion transfer 

term t± T^i^SA^+^AfiA^+lM + kc ) ^SS^P- 



two such systems described by the indices (a, b) 
we have the effective Hamiltonian for the charge sector 
of the coupled system in terms of a pseudospin ladder 
model 



j,a,b^a 



+ t± ( T i,a T lb + 



)] ,(2) 



where a,b = 1, 2 is the chain index. Having explored the 
case of U± >> V in an earlier work . we will explore 
the scenario of U± << V below. In the spin sector, cou- 
pling between the two systems leads to a S = 1/2 Heisen- 
berg ladder-type model, and has been studied extensively 
by several authors (Jills!- Interchain spin coupling turns 
out to be relevant, opening up a spin gap concomitantly 
with generating AF/dimer ordered ground states with 
spontaneously broken SU(2)/translation symmetries. 

We now analyse the charge sector. To begin, we in- 
troduce (bond-)fermion operators via a Jordan- Wigner 
transformation. In terms of these fermions, upon denot- 
ing the chains as a =|, |, we find an effective Hamiltonian 
for the ID Hubbard model with an equal-spin pairing 
term and an on-site spin-flip term 



%-,T%4 + ^X^L^j+l,a + h - c ) + n J> ( 3 ) 



3,0. 
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where t and t± are the in-chain and inter-chain hopping 
parameters respectively, \x = — It the chemical potential, 
U = —Ux the on-site (Hubbard) interaction coupling and 
V = V/4 the pairing strength. Note that while we treat 
the parameters t and V as independent parameters for 
the sake of generality, t = V = V in our original model 
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(J5J. Bosonising in the usual way (in terms of the usual 
charge p =f + J, and spin a —\ — J, variables), we obtain 
the effective low-energy bosonic Hamiltonian 

+_L /* da;Kx a ^n ff ^)) 2 + |^(a x CT (x)) 2 ] 

^7T J Act 

U f U f 

H — / dxcosf-v/S^) H — / dxcos(y/8<fin) 

zirct J 2ira J 

+ - — / dxcos(V29 p )cos(V26 a ) / dxd x (f> p (x) 

2na J 7r J 

+ — ^- [ dx cos(v / 20 (T ) cos (V28a) H / efec< 

7ra J na J 
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1 1 dxcos(V2<f> (T )cos{V20 p ) + — I dxcos(V89 p ) (4) 



where IT 



n„ 



, VpK p = v F 
-) and Va-/K a = up(l 



■^j^). Among the various cosine potentials, we have 

the usual spin-flip backscattering cos(v / 8</> -) and Umk- 
lapp cos(s/S(f> p ) terms as well as the triplet supercon- 
ducting cos(v / 20 p ) cos(v / 2# cr ) term. The chemical po- 
tential term can be absorbed by performing the shift 

<Pp —> 4>p + ^-^- x. The cosine potentials with cou- 
plings Vi, V 2 and V 3 are generated under RG by the t± 
and V terms. We find the RG equations for the various 
couplings to second-order as 
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dV 3 
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(5) 



The RG equations for the two interaction parameters 
(K p , K,j) as well as the parameter 8 = K p /j,/v p are found 
to be 



dK a 

dl 
dK, 



KlUl + V( 



dl 
dS 

u 



" - -K 2 p U 2 p J (5(l)a) 
5(l)~U 2 pMS(l)a) 



where 6(1) — 6e, a is a short-distance cut-off like the 
lattice spacing and Jo (x) , J\ (x) are Bessel functions Q . 

For repulsive interactions (U± > 0) between the bond- 
fermions, the a sector is massless and K a flows under 
RG to the fixed point value K* > 1 and 1/2 < K p < 1. 
At 1/2-filling (for the bond-fermions), the couplings U p , 
V, t±, Vi and V 2 are all relevant while U a and V3 are 
irrelevant. The competition to reach strong-coupling first 
is, however, mainly between U p , tj_ and V. We show 
below the phase diagram as derived from this analysis. 

l.l''l.21.31.4'l.5 
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K a < (1/4(1/%, + l/K„), l/4(K a + l/K a )) 
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FIG. 1: The RG phase diagram in the (K a ,K p ) plane for 
repulsive interchain interactions (U± < 0). The three regions 
K p < (1/4(1/*, + l/K p ), l/4(K a + 1 / K a ) ) , l/K a > K p > 
1/4(1/ K a + l/Kp) and K p > (l/if„ 1/4(K„ + 1/K a )) give 
the values of (K a , K p ) for which the couplings U p , t± and V 
respectively are the fastest to grow under RG. 

In the phase diagram in Fig.QJ, the three lines with 
intercepts at (K a = l,K p = 1), (K = l,K p = 0.64) and 
(K a = 1, A' p = 1/2) are the relations K a = l/K p , K p = 
1/4(1/A CT + l/K p ) and K p = l/4(K a + 1/ K a ) respec- 
tively. The regions K p < (1/4(1/ K a + 1/K p ), 1/4(A CT + 
1 / Kfj ) ) , 1 / K a > K p > 1/4(1/A CT + l/Kp) and A p > 
(1/ K a . 1 / A(K a + 1 / K a )) signify the values of A p and K a 
for which U p (in-chain Wigner charge-ordered Mott in- 
sulator), t± and V (in-chain Peierls charge-ordered Mott 
insulator of preformed bond-fermion pairs) respectively 
are the fastest to reach strong-coupling. The RG equa- 
tions for the coupling U p , the interaction parameter K p 
and the incommensuration parameter S are familiar from 
the literature on commensurate-incommensurate transi- 
tions [l^. For temperatures T » K p p,, the finite chem- 
ical potential is unable to quench the Umklapp scatter- 
ing processes, allowing for the growth of U p to strong- 
coupling. For T << K p fi, the finite chemical potential 
cuts off the RG flow of U p , freezing the Umklapp scat- 
tering processes. 

For the case of t± being the most relevant coupling, 
we find our RG equations to be a non-trivial general- 
isation of those derived in for the model of two 
coupled spinless fermion chains without the intrachain 
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V pairing term. The resulting picture then describes 
strong interchain two-particle correlations between bond- 
fermions sharing a rung. Following the analysis outlined 
in 0, 0] , we conclude that this phase is a novel insu- 
lating phase characterised by a mass gap and delocalised 
bond-fermions on rungs, resembling the orbital antifer- 
romagnetic ground state found in the spinless two-chain 
problem. This matches our finding of an orbital anti- 
ferromagnetic ground state in the strongly-coupled lad- 
der with dominant antiferromagnetic rung-couplings in 
an earlier work |l2j . Away from 1/2-filling (for the bond- 
fermions), the competition is mainly between t± and V. 
For V reaching strong-coupling ahead of t ± , the system 
is in a channel triplet-spin singlet superconducting phase 
with mobile intra-chain hole pairs; for t± reaching strong- 
coupling ahead of V , we are currently unable to describe 
in more detail the dominant instability away from the or- 
bital antiferromagnetism like insulating phase away from 
quarter-filling. 

Interestingly, while studies of ladder models have 
shown a plethora of charge and spin-ordered gapped 
phases |iLatU|, there remains the intriguing possibility 
that some of these gapped phases may themselves be sep- 
arated from one another by non-trivial gapless phases of 
finite width |l3l |. In what follows, we provide an explicit 
realisation of this scenario. The RG equations ©-© 
reveal the existence of a non-trivial fixed point (FP) for 
any value of (K p , K a ) lying the in ranges 1/2 < K p < 1, 
1 < K a < 2 + y/3 and which is perturbatively accessible 
from the trivial weak-coupling FP. We note that a sim- 
ilar non-trivial fixed point was found in a study of the 
anisotropic Heisenberg spin-1/2 chain in a magnetic field 
(lif. where the authors derived a set of RG equations 
which were very similar to those found in |l5|. Here, the 
non-trivial FP is given by 
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V* = JaK*JcK* - (K* + l fab) , ^ 




where a = 2- (K* + l/K* p )/2 ,6 = 2- (l/K* + l/K* p )/2 
and c = 2 - (K* + l/K*)/2 . Further, we can safely 
make the approximation of the renormalisations of K p 
and K„ being small at this non-trivial FP 0] . The sys- 
tem is gapless at this non-trivial FP as well as at points 
which flow to it. The trivial FP has 6 unstable direc- 
tions (Up, V, tj_, Vi, V2 and 5), 2 stable directions (U a 
and V3) and 2 marginal directions (K p and K&). The 
non-trivial FP has 5 unstable directions, 3 stable direc- 
tions and two marginal directions. Apart from the stable 
direction given by the V3 , the presence of the other two 
stable directions at the non-trivial FP indicates the ex- 
istence of a two-dimensional surface of gapless theories 
in the five-dimensional (U a , V , tj_, Vi, V2) coupling space. 
This gapless phase is the analog of the "Floating Phase" 



found in the phase diagram of the ID axial next nearest 
neighbour Ising model [l^. We present in Fig.(0| be- 
low a RG flow phase diagram which is projected onto the 
(^l) t±) plane (a similar RG flow diagram is found for the 
case of the anisotropic Heisenberg model in a magnetic 
field h in the (a, h) plane, where a is the anisotropy 
parameter) . 




FIG. 2: The RG phase diagram in the (Vi,t±) plane. The 
thick line characterises the set of points which flow to the 
intermediate fixed point (V*,tj_) shown by the filled circle. 
The thin lines show all RG flows which flow towards strong- 
coupling in the two phases I and II, characterised by the 
relevant couplings V and t± respectively. 

The regions I and II characterise all RG flows which 
do not flow to the intermediate fixed point at (V*,^). 
In region I, V± flows to strong-coupling while t± decays; 
for 1/2 < K p < 1 and K a > 1, we know from the above 
discussion that in this region, the coupling V will reach 
strong-coupling first. In region II, both t± and V\ grow 
under RG, with the coupling t± being the first to reach 
strong-coupling. Thus, the RG trajectory leading to the 
intermediate fixed point represents a gapless phase sep- 
arating the two gapped, charge-ordered phases I and II 
characterised by the relevant couplings V and t± respec- 
tively. 

For attractive interactions (U± < 0) between the bond- 
fermions, we can carry out a similar analysis. In this case, 
we can see that K p > 1 while K a < 1. Then, from the 
RG equations given above, we can see that the Umklapp 
coupling Up and V\ are irrelevant while the couplings t± , 
Ucr, V, V2 and V3 are relevant. The competition to reach 
strong-coupling first is, however, mainly between U a , tj_ 
and V. We show below the phase diagram at 1/2-filling 
for the bond-fermions as derived from this analysis. 

In the phase diagram in Fig.J2J), the three lines with 
intercepts at (K p = l,K a = 1), (K p = l,K a = 0.64) 
and (K p = l,K a = l/v3) are the- relations K a = 1/K P , 
K a = 1/4(1/ K a + 1/Kp) and K a = 1/V3 respectively. 
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1.1 1.2 1.3 1.4 1.5 

K a > (1/K P , 1/4{1/K P + 1/Ka)) 




K a < (l/i(l/K P + 1 / K& ) j 1/V3) 



FIG. 3: The RG phase diagram in the (K p ,K a ) plane for 
attractive interchain interactions (U± > 0). The three regions 
K a < (l/i(l/K a + 1/K P ), 1/V3), 1/K P > K a > l/y/3 and 
K a > (l/K p , l/4(l/K p + l/K a )) give the values of (K a , K p ) 
for which the couplings U a , t± and V respectively are the 
fastest to grow under RG. 



The regions K a < (1/4(1/K a + l/K p ),l/VS), 1/K P > 
K a > l/y/3 and K a > (1/K P , l/A(l/K p + l/K a )) signify 
the values of K p and K a for which U a (rung-dimer in- 
sulator with in-chain Wigner charge-ordering), t± and 
V (insulator with in-chain dimers and Peierls charge- 
ordering) respectively are the fastest to reach strong- 
coupling. This matches our finding of a ground state 
with in-chain Wigner charge order and rung-dimers in the 
strongly-coupled ladder with large ferromagnetic rung- 
couplings in an earlier work |l2T ]. Away from 1/2-filling 
(for the bond-fcrmions), depending on which of the three 
couplings tx, V and U a is the first to reach strong- 
coupling, the system exists either as a superconductor 
with intra-chain hole pairs (V) or a superconductor with 
rung-singlet hole pairs (U a ) or a phase reached by fol- 
lowing the dominant instability away from the orbital 
antifcrromagnetism like insulating phase (t±) but which 
we are currently unable to describe in greater detail. 

To conclude, we have studied a model of strongly corre- 
lated coupled quasi- ID systems at 1/4-filling using an ef- 
fective pseudospin ladder model with V > t, t±, U±. Us- 
ing a bosonisation analysis, we find two different types of 
charge/spin ordered ground states at 1/4-filling. Trans- 
verse bond-fermion hopping is found to stabilise a new, 
gapped (insulating) phase characterised by interchain 
two-particle coherence of a type resembling orbital an- 
tiferromagnetism 0, 0] . The spin fluctuations are de- 
scribed by a S — 1/2 Heisenberg ladder-type model for 
all cases studied here: the spin excitations are always 



massive. Away from this filling, either intra- or inter- 
chain superconductivity in a gapped spin background is 
found to be the stable ground state. We also find the 
existence of an intermediate gapless phase lying in be- 
tween two gapped, charge-ordered phases (characterised 
by the relevant couplings V and t± respectively) in the 
RG phase diagram of our model. Our analysis is espe- 
cially relevant to ladder systems like Srn- x C , a x Cu2^0 \\ 
and a — iVaV^Os or (3 — Na^.^V^O^ (a superconductor) 
which exhibit charge/spin long range order at x = and 
superconductivity beyond under pressure and/or dop- 
ingUi- 
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